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Abstract 

Copula modeling is beneficial and popular in a wide range of statistical fields, including banking, hydrology, 
geostatistics, engineering, economics, climatology, and so on. In this study, we introduce a new bivariate copula 
involving tuning parameters, polynomial and sine functions. It generalizes a previously known but understud- 
ied sine copula in the literature. We present it mathematically and graphically, and show how it can be used to 
generate new bivariate distributions with diverse support. In the remaining part, we investigate some ordering, 
stability, and dependence properties. Interesting relationships between it and the Farlie-Gumbel-Morgensten 
copula are discussed. 
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Introduction 


The notion of copula is, first of all, mathematical. A copula is a multivariate cumulative distribution function 
(CDF) whose each marginal CDF is the CDF of the unit uniform distribution. It can be used to (i) create mul- 
tivariate distribution functions with a variety of dependency structures, or (ii) study the dependence structure 
of an existing multivariate distribution independently of its marginal distributions. Reference [18] published 
the first version of copula in 1959, and the concept has received a great deal of attention since then, becoming 
the subject of numerous scientific studies. Applications of copulas to multivariate survival analysis can be found 
in Reference [9]. We may also refer specifically to Reference [8] for applications in finance, Reference [1 1] 
for applications in geostatistics, Reference [22] for applications in hydrology, Reference [20] for applications 
in engineering, and Reference [17] for applications in floods. 

As introduced in Reference [18], a bivariate copula is a bivariate function C(u, v), (u,v) € [0, 1]? that 
satisfies the following conditions: 


e Boundary conditions: for any (u, v) € [0, 1]?, 


C(u, 0)=C(0,v) =0, C(u, l)=u, C(A, =v, CA,1)=1. 


e Two-increasing condition: for any (u1, ug, v1, v9) € [0, 1]*, such that uj < ug and v1 < v9, 
C(ug, v2) — C(ug, v1) — C (u1, vg) + C (u1, v1) = 0. 


According to Reference [21], for any twice differentiable copula, this condition is equivalent to the more 
tractable condition: 8°C (u, v)/dudv > 0. 
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As a major interest in the concept of copula, the theorem of Sklar states that, if F(x, y) is a joint CDF of 
a continuous bivariate distribution with marginal CDFs given by Fy (x) and Fy(y), respectively, then there 
exists a copula C(u, v) such that F(x, y) = C(x (x), Fy(y)). Over time, a variety of copulas, or methods for 
constructing copulas, have been proposed. Examples of copula are the Ali-Mikhail-Haq copula, Joe copula, 
Clayton copula, Frank copula, Gumbel-Hougaard copula, Farlie-Gumbel-Morgensten (FGM) copula, normal- 
type copulas, cubic-type copulas, and Cuadras-Augé copula. The complete list of classical copulas can be found 
in Reference [14]. Modern developments can be found in References [10], [4], [2] and [8]. However, only a 
few copulas, to our knowledge, use trigonometric functions. We may mention Reference [1, Example 1 (point 
4)], Reference [7, Example 5], Reference [5], or the bivariate sine-type copulas developed in Reference [13]. 
In particular, the copula presented in Reference [1, Example 1 (point 4)] is defined as 


1 é 
C.(u, v) =uv +0—> sin(zu) sin(xv), (u,v) € [0, 11°, (1) 
m 


where 6 € [-1, 1] is a dependence parameter. We call it the simple PS (SPS) copula. It can modelize large 
dependences, according to Reference [1]; several measures of dependence are provided in [1, Table 1 (point 4)] 
to support this claim. Despite its undeniable interest, the SPS copula appears to have received little attention 
in the literature. The aim of this research is to highlight its characteristics and to go beyond them by proposing 
a simple extension of this copula. The proposed extension intends to add more functionalities by introducing 
two tuning parameters. It remains simple, but more flexible than the SPS copula in terms of modulation 
of the trigonometric functions, and better suited for the development of new statistical models in a variety of 
fields. We illustrate this claim by studying some of its properties, such as three equivalent expressions involving 
different functions, the generation of new interesting bivariate distributions, a possible multivariate extension, 
some ordering properties, including one involving the FGM copula, a stability property, the determination of 
various measures of dependence, and diverse tail dependence characteristics. 

All these aspects are examined in detail, with the following organisation. The proposed copula is defined in 
Section . Section focuses on its major properties. Finally, in Section , the research findings are summarized. 


New polynomial-sine copula 


Presentation 
The following proposition presents a new simple copula derived from Equation (1). It is based on tuning 


parameters, polynomial and sine functions. 


Proposition 1 Let (m,n) € [1,+00)?, 0 € [-1, 1], and C,(u,v), (u,v) € [0,1]? be the following bivariate 
Junction: 


C,(u,v) = uv + j [sin(ru)]” [sin(zv) |”. (2) 
n?mn 
Then C,(u, v) is a copula. 


Proof. Let us check that the boundary and two-increasing conditions defining a copula are satisfied. 


e Boundary conditions: since sin(0) = 0, it is immediate that C, (u, 0) = C,(0, v) = 0. Now, since sin(x) = 
0, we get the other conditions: C. (u, 1) =u, C,(1, v) =v, and C,(1, 1) = 1. 
e Two-increasing condition: since C, (u, v) is twice differentiable and cos(x), sin(x), 6 € [-1, 1], itis clear 
that 
3?C,(u, v) 
Oudv 


The two-increasing condition is established. 


= 1 + 8 cos(ru) cos(av)[sin(zu)]""![sin(av)]""! > 1 - || > 0. 
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This ends the proof of Proposition 1. o 





To our knowledge, the copula defined as Equation (2) is new in the literature. It is reduced to the SPS 
copula for m = n = 1. We call it the polynomial-sine (PS) copula. Note that, for @ = 0, we have C, (u, v) = uv; 
C, (u, v) is thus reduced to the “independent” copula. 


The PS copula density is given as 


0?C, (u, v) 
Oudv 


Cx (U,V) = = 1+ 8 cos(ru) cos(av)[sin(xu) |"! [sin(av)]""1. 


It is clear that m and n affect the possible shapes of this function in a sophisticated manner. Basically, for 
m € (1, +00) orn € (1, +00), we have c,(0, 0) = c.(1, 1) = 1, which is not the case when m = n = 1, where 
c,(0, 0) =¢,(1, 1) = 14+ 6, which strongly depends on 8. 


For m = n, C,(u, v) = C.(v, u) for any (u, v) € [0, 1]?; the PS copula is symmetric. It is not symmetric for 
m + n. For @ € [0, 1], we have C,(u, v) > uv, implying that any random variables X and Y having C,(w, v) 
as copula are positively quadrant dependent (PQD). For 6 € [-1, 0], the reversed inequality holds; they are 
negatively quadrant dependent (NQD). The reflected copula associated to C,(u, v) is given by C.(u, v) =ut 
v-1+C,(1-u, 1-0), so 


C.(u,v) =utv-1+(l-u)(l-v)+ e- iodi — u))|” [sin(a (1 — v))]” 


C 
n2mn 


[sin (xu)]” [sin(zv)]” = C(u, v). 





=uv+0 5 
n mn 


This shows that the PS copula is radially symmetric. As an additional property, like any copula, the PS copula 
satisfies the Fréchet-Hoeffding inequalities: for any (u, v) € [0, 1]?, we have max(u +v — 1, 0) < C, (u, v) < 
min(u, v). 

More secondary, from the analytical standpoint, we can express C, (u, v) in different manners. Two of them 
are presented below. 


e First, by introducing the normalized sinc function defined by sinc(x) = sin(ax)/(zx), we have 


c 
g+” 2 


C. (u, v) = uv {1 + 0 — u" py" [sinc(u)]” sine(]"| : 
mn 


e Second, the following elementary trigonometric formula: sin(zw) = 2 sin((2/2)u) sin((7/2)(1 - u)), 
gives 


C,(u, v) = w+ 0m [sin (Fu)]” [sin (5a -u))|" [sin (Fo) [sin (Fc -v))]". (3) 


Some similarities between this representation in its simplest form and the FGM copula will be discussed 
later. 


We complete this description with the three-dimensional and contour plots of the PS copula for selected values 
of the parameters 0, m and n. 


3D plot of the PS copula with @=1, m=1 and n=1 








0.0 1.0 


Figure 1: Three-dimensional and contour plots of the PS copula defined with 6 = 1 and m =n = 1. 


3D plot of the PS copula with 0@=-1, m=1 and n=1 
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Figure 2: Three-dimensional and contour plots of the PS copula defined with 6 = -1 and m =n = 1. 


3D plot of the PS copula with @=1, m=30 and n=1 
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Figure 3: Three-dimensional and contour plots of the PS copula defined with 6 = 1, m = 30 and n = 1. 
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Contour plot of the PS copula with 0=1, m=1 and n=1 
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Contour plot of the PS copula with 0=-1, m=1 and n=1 
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Contour plot of the PS copula with @=1, m=30 and n=1 
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Figure 4: Three-dimensional and contour plots of the PS copula defined with 6 = —1, m = 30 and n = 1. 


From Figures 1, 2, 3 and 4, we see that the PS copula has the aspect of a triangle that deforms itself 
depending on the values of the parameters. Also, a more or less pronounced bump is sometimes observed, 


certainly due to the effect of the sine terms combined with the power parameters m and n. 
Remark 2 As new facts on the existing SPS copula, it can be expressed in terms of cosine functions as 
1 
C,(u, v) = uv + 0-3 [cos(au) cos(zv) — cos(a (u +v))], 
m 
or in terms of the normalized sinc function as 
C, (u, v) = uv |1 + 8 sinc(u) sinc(v)] , 


or in terms of more sine functions as 





4 
C, (u,v) = uv + 0-3 sin (5x) sin (5v) sin (Fa — u)) sin (Fa -v)) ; 
We mention that, as the parental PS copula, it covers the independent case, satisfies the PQD property for 0 € [0, 1] and 
the NQD property for 0 € [-1, 0], and ts radially symmetric. Furthermore, it is symmetric. It was previously plotted in 
Figures 1 and 2 with 0 = 1 and 0 = -1, respectively. Last but not least, the SPS copula density is quite simple; it is given 
as 


Cx(u, v) = 1+ 0 cos(zu) cos(zv). 


Some generalized bivariate distributions 


On the other hand, the PS copula can be utilized to create generalized bivariate distributions. A few examples 


are given below. 


e The PS normal (PSN) distribution: Let u € R, oœ € (0, +œ) and ®(z; u, o) be the CDF of the nor- 
mal distribution with mean y and standard deviation ø, i.e., B(x; u, 7) = (1/(o V27)) (ie exp(-(t - 
u)? /(207))dt, x € R. Then, the following CDF defines the PSN distribution: F (x, y) = C,(®(2; u1, 01), (y; u2, w2)), 
(x,y) € R?, with (u1, u9) € R? and (a1, oy) € (0, +00)?. In an expanded form, we have 


1 . m $ n 
F(x, y) = ®(a; p1, 01) PY; wo, o2) + OE [sin (7®(x; 41, @1))]” [sin (7B; u2, o2))]" , 


(x,y) € R?. 
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As an example of representation, Figure 5 shows the three-dimensional and contour plots for the PSN 








distribution defined with the following arbitrary configuration: 6 = 0.5, m = 8,n = 1, wy = u2 = 0 and 


Oo] =o9g=1. 


Example of 3D plot of the PSN distribution Example of contour plot of the PSN distribution 

















Figure 5: Three-dimensional and contour plots of the PSN distribution defined with 6 = 0.5, m = 3, n = 1, 
My = H2 =O and oj = 09 = 1. 


e The PS Weibull (PSW) distribution: Let G(x; 8,4) be the CDF of the Weibull distribution, i.e., 
G(a;B,A) = 1 — exp(-Ax*) for x € (0, +œ), and G(x; B,A) = 0 elsewhere. Then, the follow- 
ing CDF defines the PSW distribution: F(a, y) = C.(G(x; B1, 41), G(y; B2, A9)), (x,y) € RÊ, with 
(B1, B2, A1, 42) € (0, +00)*. In an expanded form, this gives 


F(x, y) = (1 = exp(-ai24")) (1 - exp(—agy**)) 


+ I [sin (x exp(—A To [sin (x exp(-Aay?*))|" , (x,y) € (0, +00)”, 
and F(x, y) = 0 elsewhere. 


Figure 6 displays the three-dimensional and contour plots for the PSW distribution defined with the 





following arbitrary configuration: @ = 0.5, m = 3, n = 1, a, = 0.5, Ag = 1.5, and £1 = Bo = 2. 


Example of 3D plot of the PSW distribution Example of contour plot of the PSW distribution 























Figure 6: Three-dimensional and contour plots of the PSW distribution defined with @ = 0.5, m = 8, n = 1, 
A, =0.5, Ag = 1.5, and Pi=Po= 9. 
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e The PS power Lomax (PSPL) distribution: Let H (x; «œ, 8, A) be the CDF of the power Lomax dis- 
tribution, i.e., H(z;a, B,A) = 1 -—A%(A + 28)-* for x e (0,+00), and H(a;a, B,A) = 0 else- 
where. See Reference [15]. Then, the following CDF defines the PSPL distribution: F(x,y) = 
C.(H (a; 01, B1, 41), H (Y; 09, B2, A9)), (x,y) € R?, with (a1, a2, B1, B2, 41,42) € (0,+00)°. In 


an expanded form, we have 


F(z,y) = (1 aah; +281) (1 sA a + yb2)-aa) 
+ g- [sin (nay (Aq +a) |" [sin (raa +y) |" , (x,y) € (0, +00)”, 


nmn 


and F(x, y) = 0 elsewhere. 


Figure 7 displays the three-dimensional and contour plots for the PSPL distribution defined with the 





following arbitrary configuration: 0 = -1, m =n = 1l, a1 = 1 = 2, aq = Bg = 3, anda; =Ag= 1. 


Example of 3D plot of the PSPL distribution Example of contour plot of the PSPL distribution 





















































Figure 7: Three-dimensional and contour plots of the PSPL distribution defined with 0 = -1, m =n = 1, 
Qa) = Pi =2,a9 = Bg = 3, and 41 = 9 = l. 





Other examples can be presented in a similar manner. To our knowledge, these three bivariate distributions 
have not been examined in a statistical modelling scenario, and have a certain interest in this regard. This, 
however, demands further investigations that we leave for future work. 


Multivariate case 


We end this part by a corollary about a natural multivariate extension of the PS copula. 


Corollary 1 Let d be a strictly positive integer, mı, mo, ..., mq € [1, +0)4, 0 € [-1, 1], and C, (u1, u9, . . . , ua), 
(u1, u9, ... , uq) € [0, 1]? be the following multivariate function: 





d d 
1 
C. (u1, Ug, ..., Ud) = | | +o | | [sin (rue) ]"*. 
k=1 mt [Ig Mk pat 


Then C,(u1, ug, ..., uq) is a copula. 


The proof can be conducted as the proof of Proposition 1. 
The multivariate version of the PS copula can be called the multivariate PS (MPS) copula. One can use it 
to create generalized multivariate distributions, and multivariate statistical models as well. 
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Specific properties 


The PS and SPS copulas satisfy a number of special properties, which are detailed in this section. 


Ordering properties 


Some ordering properties satisfied by the PS copula are now being examined. The following result is about a 
direct hierarchical order involving the PS and SPS copulas. 


Proposition 3 Let C,(u, v) = C. (u, v; m, n) be the PS copula with the specification of m and n, and C,(u, v; 1, 1) be 
the SPS copula. For any 6 € [0, 1] and any (u,v) € [0, 1]?, we have 


C,.(u, v;m,n) < C(u,v; 1, 1). 
For any 0 € [-1, 0], the reversed inequality holds. 


Proof. For any m € [1, +00), we have [sin(au)]”~! < 1 < m, implying that [sin(au)]”/m < sin(zu). There- 


fore, for 6 € [0, 1], since all the involved quantities are positive, we have 





1 
C(u, v;m, n) =uv +0 [sin(xu)]” [sin(zv)]” < uv + 0 —z sin(xu) sin(zv) = C,(u, v; l, 1). 
m 


7 
nmn 


The desired result is proved. The inequality is reversed by considering a negative value of 0. This ends the 





proof of Proposition 3. o 


The following result shows a double inequality involving the PS copula and a function connected to the 
FGM copula, and some of its generalizations. 


Proposition 4 With the notations used for the PS copula, let us set 
A(u, v; 0) = uv + Ou™(1—u)™v"(1—v)", (u,v) € [0, 1]?. 


For any 6 € [0, 1] and any (u, v) € [0, 1]?, the PS copula satisfies the following double inequality: 


m+n 





Aine 


grz(mtn-l) 
| Clu») < A (n056 ) 


n?mn Qm mn 


For 0 € [-1, 0], the reversed double inequality holds. 


Proof. Let us have in mind the representation of the PS copula in Equation (8). The following double in- 
equality holds: for y € [0, 7/2], (2/7)y < sin(y) < y. Therefore 


m m_n n . [T Mmf, (7 al EEE 1s nT [T n 
u” (1 —u)”v” (1 -— v)" < [sin (Zu)] [sin (Fa - »))| [sin (5v) [sin (Sa -v))| 
gen) P in ae i 
S a” (1-u)”"v"(1 - v)”. 
We thus obtain the desired result by applying this double inequality to Equation (8). The double inequality is 


reversed by considering a negative value of 6. This ends the proof of Proposition 4. Oo 


Let us now discuss some ordering relationships between the PS and FGM copulas. First, the FGM copula 
is defined by 
Co(u,v) = uv + Ouv(1—u)(1—v), (u,v) € [0, 1]?, 


with @ e [-1, 1]. Hence, for m = n = 1, we have A(u,v; 6) = Co(u, v). Therefore, for some values of 6, m 
and n, Proposition 4 provides direct ordering properties between the SPS and FGM copulas. In addition, for 
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any (m, n) € [1, +00)”, we have A(u, v; 0) < Cs (u, v). Therefore, for some values of 8, m and n, Proposition 4 
shows that the PS copula can be compared to the FGM copula. In addition, ordering properties between the 
parental PS copula and some general forms of the FGM copula involving powers of u, 1 — u, v and 1 — v can 
be established, such as the one proposed in Reference [12] by putting the power parameters equal. 


Stability property 


In this subsection, we investigate a special stability property satisfied by the PS copula. It involves a special 
product presented below. Let C\ (u, v) and C? (u, v) be two copulas. Then we define the copula product of 
C (u, v) and C® (u, v) by 





, 1 acO) (2) 
(CO 4C®)(u, v) ad ge eee ey, (4) 
0 Ov Ou 


It is proved in Reference [14] that (C * C))(u, v) is also a copula. In the following result, we establish a 
stability property of the PS copula with respect to the product above. 


Proposition 5 For any q € (0, +00), let us set 
1 
(q) = | [sin(ax)]%dx. (5) 
0 
Let C,(u, v) = C,(u, v; 0) be the PS copula with the specification of 0, ce (u,v) = C(u,v; 01) and cP (u, v) = 
C, (u, v; 09) with 01, 09 € [-1, 1]?. Then (Cc «CY (u, v) = C, (u, v; 0x), where 
0x = 0102 [I(m +n — 2) —- (m + n)] . 

Thus, the copula product of two PS copulas is also a PS copula. 


Proof. Based on Equation (1), we have 


ðC. (u, v) 


=u+t oL cos(ru)[sin(zu)]”} [sin(srv) ]” 
Ou nn 


and 
aC (u,v) 


1 
= u + 0— cos(xv)[sin(xu)][sin(v)]""!. 
ðv nm 
Therefore, based on the copula product in Equation (4), we have 


(CP x) (u,v) 


i 
0 


= uv + 092 


l cos(t)[sin(zt) |”! [sin(av)]" | dt 
nn 


v+ 09 





u +01 L cos(t)[sin(zu)]”[sin(zt) |"! 














mn 


1 
alsin)" f cos(7t)[sin(at)]” tdt 
0 


1 
+4) olsin)" f cos(xt)[sin(st)]""!dt 
mm 0 


+6169 





1 
[sin(ru)]" [sin(7v)]” I [cos(mt)]? [sin(x1)]"*"-2di 


G 
m2mn 


=uv+0 +0 +0102 





5 [sin(au)]” [sin(arv)]” [S(m +n — 2) - I(m + n)]. 
n?mn 


We recognize the PS copula with dependence parameter 6, = 0109 [J(m +n — 2) —- (m +n)], ending the 
proof of Proposition 5. o 
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In particular, if CO (u, v) and cP (u, v) are two SPS copulas with parameters 6; and 69, then (Cc * 
C)(u, v) is the SPS copula with parameter 0182/2. A similar stability property also holds for the FGM 
copula (see Reference [19]). We end this part by discussing the integral 3(q) presented in Equation (5). Even 
if there is not a simple expression for this integral, it can be calculated for any reasonable value q via standard 
integral or numerical techniques. In order to illustrate this claim, Table | presents some of its values for selected 
integer values of q. 


Table 1: Values of 3(q) for several integer values of q 



























































3(q) | 1 2/1) 4/3] 16] 5 | 82 | 35 | 256 | 63 
i m | 2 | Ən |8| 15r | 16 | 35m | 128 | 315 | 256 
q 11 12 13 14 15 16 17 18 19 20 





512 | 231 | 2048 | 429 | 4096 | 6435 65536 | 12155 | 181072 | 46189 
6987 | 1024 | 80087 | 2048 | 6435x | 32768 | 1098957 | 65536 | 230945r | 262144 





3(q) 









































As basic properties, it is decreasing with respect to q. Note that 3(q) will also be involved in the dependence 
properties of the PS copula. 


Dependency properties 


In the bivariate scenario, measures of dependence are frequently used to describe a complex dependence struc- 
ture. This part is devoted to some measures of dependence related to the PS copula. 

Let X and Y be continuous random variables whose copula is denoted by C(u, v). Then the dependence 
of X and Y can be evaluated by several measures, such as the Kendall or Spearman measures. The Kendall 
measure is defined by 





1 1 
=j -4 f | OC(u, v) OCU) 5 ae (6) 
0 0 Ou Ov 


and the Spearman measures is specified by 


1 1 
p= 12 f Fi C(u, v)dudv — 8. (7) 
0 0 


For alternative expressions and details on these measures, see Reference [16]. 


When the PS copula is considered, these two measures are determined in the next result. 


Proposition 6 The Kendall and Spearman measures of dependence related to the PS copula are given by 


T= 884 





z S(m)S(n),  p=120— 
n?mn n?mn 





(m) I (n). 


10 
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Proof. Based on Equation (2), we have 


OC, (u, v) 
Ou 





=v+0 L cos(zu)[sin(zu) |”! [sin(srv) ]” 
nn 


and 


OC,(u, v) 
ðv 


Therefore, in view of using Equation (6), we calculate 


ôC, & v) OC,(u, v) 
a f a0 dudv 
i [ A i m= [si ‘| l : mps B 
= v + 0 — cos(ru)[sin(ru)]” [sin(zv)]” | (u + @— cos(rv)[sin(ru)]” [sin(av)]"* | dudv 
o Jo nn =m 
1 1 1 
7 y udu! ane +0 f vcos(7v) [sin(ra) "1d f sin(a) Pd 


1 1 
o+ p u cos(ru)[sin(m)]"= du f [sin(o)] "do 


1 1 
+ ee f cosu) [sin(zu)]?" "da i cos(7v) [sin(ao)] =de] ‘ 


nmn 


pA 





=u+60 L cos(zv)[sin(zw)]” [sin(zv) 
nm 











By an integration by part, it comes 


1 
1 v cos(zv) [sin(zv)]""!dv = s a: 
0 nn 


The other integral terms can be calculated directly; we obtain 


ie l a 
0 Jo ðv "=I mn 














i 


Thus, by Equation (6), the Kendall measure is obtained as 


| =80 


For the Spearman measure as described in Equation (7), the calculus is more direct; it follows from Equa- 


tion (2) that 


1 1 1 1 
a i Cou, oldu =| fo au] | I odol +0 | f fsin(ru) da | f [sin(ro)]"d 








1 
ra1-4(7 








1 1 
=—+ O amn 2 
4 
Therefore i 
= 12|- 3= 
P | 7 +0 |- 
This ends the proof of Proposition 6. o 





By taking m = n = 1 in Proposition 6, we rediscover the result for the SPS copula described in [1, Table 1 
(point 4)]. That is 
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It can be noted that 32/74 = 0.8285114, and 48/74 ~ 0.4927671, which give the following approxima- 
tive dependence domains: t € [—0.338, 0.388] and p € [-0.5, 0.5]. These domains extend the dependence 
domains associated with the FGM copula which are r € [-0.23, 0.28] and p € [-0.34, 0.84]. Thus, the 
FGM copula does not allow the modeling of high dependences, whereas the SPS copula does. The SPS cop- 
ula is preferable in this aspect. It is also preferable to some modern extensions of the FGM copula, as those 


presented in References [10] and [8]. The dependence domains of the more general PS copula are: 


1 
rm) 3(n), 12—— S(m)3(n) i 
n?mn n?mn 








TE -8——3(m)3(n), 8 Sma) i pE |-12 
n mn 


n?mn 

Since 3(m) and 3(n) decrease as m and n increase, the maximal dependence domains are obtained for 
m =n = 1. The addition of m and n to the SPS copula provides more flexibility in the possible shapes of the 
main functions, but does not increase the range of the dependence domains. 


Tail dependence 


The level of dependency in the upper-right quadrant tail and lower-left quadrant tail of a bivariate distribution 
is referred to as tail dependence. It is a concept that is relevant to the research of extreme value reliance. Let 
X and Y be continuous random variables. Then the tail dependency between X and Y turns out to be a copula 
property; the amount of tail dependence is invariant under strictly increasing transformations of X and Y. See 
Reference [6]. By denoting C(u, v) the copula of X and Y, the following notions are adopted: 


¢ Let us consider 
> 1L-2u+C(u, u) 
au = lim ———————_. 
ul l-u 
If Ay exists, then C(u, v) has an upper tail dependence if Ay € (0, 1], and no upper tail dependence if 
Ay = 0. 


e Let us consider 


C 
irsin n 
u0 u 


If Az exists, then C(u,v) has a lower tail dependence if Az € (0, 1], and no lower tail dependence if 
Ay = 0. 


Extreme value theory makes considerable use of these two measurements. Here, we investigate them for the 
PS copula in the next result. 


Proposition 7 The measures of tail dependence related to the PS copula are given by Ay =A, = 0. 


Proof. Owing to Equation (2), when u — 1, the following equivalence holds: 

















1-2u+ C.(u, u) =1-u+0 1 [sin(au)]”"*" = al = u) l+ gl gmn-2(q = wor 
l-u n?mn l-u mn 
Hence ay = 0. 
Moreover, when u — 0, we have 
Cx (u, u) =u+0 l [sin(au) ]"*" sull gL prtn-2 ymin- 
u nmn u mn ` 
This entails that A; = 0. The proof of Proposition 7 is now complete. o 


Since Ay = A, = 0, the PS copula has no upper tail dependence or lower tail dependence. 
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Summary with discussion 


In this article, we have provided a new bivariate copula extending the one in Reference [1, Example 1 (point 
4)]. It is based on tuning parameters, polynomial and sine functions. We call it the polynomial-sine copula. We 
have described it in detail, and illustrated how it may be used to produce new bivariate distributions with various 
supports. In addition, we have discussed some connections with the Farlie-Gumbel-Morgensten copula, as well 
as various stability and dependence properties. Among its main features, the polynomial-sine copula has the 
following merits: (i) it has a simple mathematical structure, (ii) it is connected, in some senses, to the Farlie- 
Gumbel-Morgensten copula and some of its extensions, (iii) it may have a dependence domain wider than the 
one of the Farlie-Gumbel-Morgensten copula, and other well referenced copulas, and (iv) it has no upper tail 
dependence or lower tail dependence. Hence, the proposed copula is believed to be beneficial in modeling real 
data sets by practitioners. Meanwhile, the author is currently employing the polynomial-sine copula on some 
real data sets, with the findings to be published in a future study. 
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